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Abstract 

\q . We report detailed results on ionization in metallic quantum-dot (QD) nanorings described by 

the extended Hubbard model at half filling obtained by exact numerical diagonalization. In spite 



o 



this scarcity to a hidden quasi-symmetry, generalizing thereby similar results on optical absorption 



of very strong electron correlations, the ionization spectra are astonishingly scarce. We attribute 
U 
CO 

recently reported [I. Baldea and L. S. Cederbaum, Phys. Rev. B 75, 125323 (2007); 77, 165339 
(2008)]. Numerical results indicate that this hidden quasi-symmetry of the extended Hubbard 
model does not evolve into a true (hidden) symmetry but remains a quasi-symmetry in the case 
of the restricted Hubbard model as well. Based on the observation on the number of significant 

ionization signals per each spatial symmetry, we claim the existence of a one-to-one map between 

> ■ 

q^ ' the relevant ionization signals of the correlated half-filled nanorings and the one-hole and two- 

VO ■ hole-one-particle processes possible in the noninteracting case. Similar to the case of optical 

r^ ', absorption, numerous avoided crossings (anticrossings) are present in the ionization spectra, which 

O 

often involve more than two states. The present results demonstrate that ionization could be a 



useful tool to study electron correlations in metallic QD-nanoarrays, providing information that is 



complementary to optical absorption. 
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I. INTRODUCTION 



Metallic quantum dots (QDs) can be fabricated and assembled in extended regular arrays 
by means of modern nanotechnologies [l|, y, y|, |j, la, la, LZ|- The salient feature of such 
artificial nanostructures is the fact that their electronic properties can be smoothly tuned 
in wide ranges by factors that can be easily controlled experimentally: dot diameter 2R, 
interdot spacing D, and gate potentials. It is this tunability of nanostructures that makes 
the most important difference from ordinary molecules and solids. In isolated QDs with 
sizes of a few nanometers, typical for the nanograins of silver prepared by Heath' group 



q, a a q, a a 



101 ]. the single-electron levels, characteristic for quantum boxes or 



"artificial" atoms Jjj , are well separated by energies ~ 1 eV. A small fraction of "valence" 
electrons, which occupy the highest "atomic" orbitals, becomes delocalized over the whole 
nanostructure, if the QDs are sufficiently densely packed and the adjacent electronic wave 
functions overlap. 

For nearly touching Ag-QDs (d = D/2R > 1), correlations of these delocalized electrons 



are negligible: 
approximation 



;he single-particle, molecular-orbital (MO) picture represents a very good 



12 



13 



14j . This description progressively worsens for more and more distant 



QDs, and completely breaks down towards the end of the <i-range of experimental interest 



(d S, 2). Electron correlations become more and more important beyond d > 



14 



1.4 |l2 



m 



151 ] . The attractive fact to study such QD-nanostructures is that by varying d one can 



drive continuously a many-electron system between weak and strong correlation regimes. 
Obviously, this route from weak to strong correlation regimes is inconceivable in ordinary 
molecules and solids. 

The present paper is devoted to metallic QD nanorings. To our knowledge, rings of 
metallic QDs have not yet been assembled so far. For this purpose, the use of a mask in the 
preparation of a two-dimensional nanoarrays l|, y, |j, |5|, |7|, |8|, |9|, llOj is conceivable. As shown 
recently 13j, ll4j , electron correlations in these metallic QD-nanorings turned out to possess 
an intriguing character. In spite of very strong correlations, the optical absorption spectra 
are astonishingly scare (almost monochromatic for closed shells 13|) and can be rationalized 



within a single-particle description. It would be interesting to further investigate to w 
extent the hidden quasi-symmetry suggested in the context of optical absorption [13j, 



hat 



ii 



manifests itself in other properties. Along this line, in the present investigation we shall focus 



our attention to (photo)ionization, a pro p ert y for which electron correlations are known to 
be important lg, [l7|, [l8|, [l9|, 120, |2l|, |22j, [23|, |24| . The important role played by electron 



correlation in the ionization of QD nanorings became evident in the preliminary study of 



Ref. 



12 



12 



In this sense, the present work in an extension of the investigation of Ref 
where only the lowest energy process was considered. This is the so-called HOMO process, 
because within the single-particle picture it would correspond to removing an electron from 
the highest occupied molecular orbital. The fact demonstrated in Ref. [l2j, that the molecular 
orbital picture of ionization completely breaks down in QD-nanorings even for the HOMO 



processes 



see Refs. 



is very interesting^ 



mU,U,U, 



2d, 



21 



because it is contrary to ordinary molecules, where it holds; 
L|. This represents a further motivation for studying ionization 
in QD-nanorings. 

To anticipate, the results for ionization we are going to present below fully support the 
existence of a hidden quasi-symmetry in the strongly correlated metallic QD-nanorings 
described within the extended Hubbard model found in optical absorption [13, llJ]. Still, 
because we cannot present an analytical proof on the existence of a hidden quasi-symmetry, 
we believe that it is important and beneficial for further investigations to amply document it 
by detailed exact results obtained by means of extended numerical calculations. Because the 
numerical study is rather intricate, in Sec. [TT1 where the model of QD-nanorings is exposed, 
we also provide significant details, which were not given in the earlier studies. In the next 
Sec. IIHI we analyze the electron correlations from a standpoint different from that adopted 



in Refs 



.0 



13 



14j . e. g., by considering the quantum entanglement. In Sec. IIVI the results 



on ionization spectra will be presented. We shall examine in detail the cases of nanorings 
consisting of six and ten QDs separately, in Sees. IIVAI and IIVB| respectively. The two 
main effects of electron correlations, the avoided crossings and the scarcity of the ionization 
spectra will be considered in Sees. [V] and IVlj respectively. Conclusions will be presented in 
the final Sec. ETTl 

II. MODEL AND METHOD 



To describe the valence e' 
series of previous studies [l2 



ectrons of metallic QD-nanorings, we shall utilize, similar to a 



jiliilia 



251 ] . an extended Hubbard (or Pariser-Parr-Pople) 



model Hamiltonian 



N 



H 



h,ff a i+\,o + a 



l+l,cr a l,a 



1=1 <r=U 

N 



(1) 



1=1 



,t 



where, a (a') denote creation (annihilation) operators for electrons, ni y(T = a' la .a la , ni 
n i,] + n>i,i) The ideal situation assumed in Eq. (pQ), where the model parameters are site 
independent, can be considered a reasonable approximation in view of the narrow size dis- 
tributions (~ 2 — 5%) in the arrays of Ag QDs assembled by Heath's group [l|, |2j, |4j, la, LZ|- 



As previously discussed [12 



131 ]. such a weak disorder does not significantly alter the results 



obtained by assuming site-independent parameters. Importantly for subsequent considera- 
tions, the point symmetry group associated with the model ([1]) of nanoring is Djvti- 

^he model parameter entering Eq. ([1]) has been analyzed earlier in literature [2j, |7|, ll2|, 
271 . l28l . |29J and will be therefore no more repeated here. Importantly, the fact that 



26 



the interdot separation d = D/(2R) (measured between QD centers) can be continuously 
varied in the range 1. 



parameter tuning [1 



._0 < d < 1.85 by means of a Langmuir technique allows a wide 

nnfTQ n 

, |2|, U, [5|, 17J • For concrete values, see Fig. 1 of Ref . [14|. At the ends of the 
aforementioned <i-range, the ratio of the on-site repulsion U to the free electron bandwidth 
4to varies from values U/(4t ) <C 1 to f//(4t ) 3> 1. Therefore, a crossover from a weakly 
correlated system to a strongly correlated one by varying d can be reached. 

Ionization is usually described as ejection of an electron from a "molecular" orbital (MO). 
Therefore, besides electron operators ai jff for "atomic" orbitals, it is also useful to introduce 
operators of MOs (or Bloch states), Ck, a = 1/yiV ^ aj )0 - exp(— 27rkli/N). An ionization 
signal is characterized by an ionization potential Ek,i = (i$!k,i\H\^k,i) ~ (*&\H\®) anc l a 
spectroscopic factor w 

w K i = \{^ Ki \c Ka \<$>)\ 2 . (2) 

Here, k denotes the wave number of the Bloch state (MO) out of which the electron is 
removed and ^^ (i = 1,2,..., see below for notation) are eigenstates of the ionized nanor- 
ing. We shall only consider ionization at zero temperature, that is, the nonionized system 
is initially prepared in its ground state $. The quantities w^i are important, since the 



photoionization signals are proportional to them 22l.l23l. |24|. 



Symmetry considerations play an important part for ionization in nanorings, similar to 



optical absorption 13j, Il4j . In the latter case, the counterpart of the matrix elements in 
the r.h.s. of Eq. (j2]) entering the coefficient of optical absorption are the elements of the 
dipole operator P. For the case of closed shells considered here, we have always found that 
the ground state \l/ possesses A\ g symmetry. In this case, because the dipole operator P 



possesses E\ u symmetry [13|, ll4j, only i?i u -excited eigenstates can be targeted in optical 



absorption. To target eigenstates with symmetries different from E\ u by zero-temperature 
optical absorption, one needs to consider open shell rings, characterized by ground state 
symmetries T$ ^ A\ gi as done in Ref. [14|. The situation is different for ionization. The 



MOs span many irreducible representations Tj, of the point group Djv/i |30j. Eigenstates 
^fe,i; ^k,2, ^k,3, • • • of all MO-symmetries (specified by the value of k) can be targeted by 
ionization even with an initial neutral Ai s -state $ (because 1^ ® Tk D A\ g ). To this, it 
suffices to consider nanorings with closed shells. For the sake of simplicity, we shall restrict 
ourselves here to the prototypical case of closed shells, the half-filling case. 

Numerical results will be presented for half-filled nanorings with six and ten QDs. For 
the former case, all eigenstates of Hamiltonian (TjQ) can be computed exactly by numerical 
diagonalization: the dimension of the total Hilbert space is equal to 924 for the neutral ring 
and 792 for the ionized ring. The neutral ground state is a singlet, and one can carry out 
calculations in subspaces with total spin projection S z = and S z = +1/2, respectively. The 
corresponding dimensions are then reduced to 400 and 300, respectively. The fact that all 
eigenstates of Hamiltonian (JTJ can be computed exactly enables us a very detailed analysis of 
all relevant physical aspects. For ten QD-nanorings, the total dimension is equal to 184756 
in neutral rings and 167960 in ionized rings. Similar to the above case, the restriction to 
the subspaces with S z = and S z = +1/2 is also possible, which results in the reduced 
dimensions of 63504 for neutral rings and 52920 for ionized rings. For such dimensions, 
the Lanczos algorithm has been applied. By running three times the Lanczos procedure we 



are able to compute the ionization energies and the spectroscopic factors individually |31 | . 
The difference from calculations for optical absorption [13j, [14], |32j, |33|] is that the first two 



Lanczos runs are carried out for the nonionized ring, while in the third run one must choose 
the normalized vector Cfc )CT |$) as starting Lanczos vector separately for each symmetry T^ 
and perform calculations for the ionized ring. 

It is worth noting at this point that the continued fraction algorithm, familiar in con- 



densed matter physics [34], |35j, |36|], does not suffice for the present purpose, since it only 
allows calculations of convoluted spectra. The phenomenon of avoided crossing, often en- 
countered in the study on optical absorption in QD nanorings [14J and also important for 
ionization (see Secs. HVA~lHVB[ and|V|), represents a typical situation where the information 
on convoluted spectra obtained by this algorithm is totally insufficient. 

Especially for the case of larger nanorings, for which the Lanczos algorithm has to be 
applied, the following sum rule, straightforwardly deduced from Eq. ([2]), 

y^Wfc.i = ($|4 )<T c fej(r |$) = n k>(T . (3) 

i 

turns out to be very useful as check of numerical calculations. 

At the end of this section, we note that, similar to earlier works [13j, [l4j, we safely rule 
out any spurious contribution, e. g., those of higher-spin eigenstates. Because the ket state 
Cfc,jJ$) entering Eq. ([2]) is an exact eigenstate of both the total spin projection and the total 
spin, S z = S = 1/2 [37|, only the spin doublet eigenstates VP^j can possess nonvanishing 
spectroscopic factors w k ,i- For the ionized six-QD nanorings, out of the total 300 eigenstates 
with S z = 1/2 we only retained the 210 eigenstates with S z = S = 1/2. In the calculations for 
the ionized ten-QD nanorings, the starting Lanczos vector possesses the correct spin (see 
above), but accumulated numerical errors during the Lanczos iterations could yield Ritz 
vectors containing unphysical admixtures with other eigenstates. We can unambiguously 
exclude this possibility: we have stored all the relevant Ritz vectors (usually at most 100 
vectors of dimensions 52920) and checked that they are very accurate eigenvectors with 
correct spin, and exhaust Eq. (3). Their dispersion was always found to be very small, 
comparable to that of the neutral eigenstate $. In particular, it is at least 3-4 orders of 
magnitude smaller than the energy difference at avoiding crossings (see Sec. M) . 

III. ELECTRON CORRELATIONS 

The fact that electron correlations are very important for Ag-QD nanorings was demon- 



strated in a series of previous works. Details can be found in Refs. 



12 



m, 



14 



15 



and will be 



therefore largely omitted here. One can inspect, e. g., the MO occupancies; see Fig. lb of 
and Fig. 2 of Ref. [14| for ten QD-nanorings. At half filling, up to say, d < 1.4, one 



Ref. 



13 



can clearly distinguish lower, occupied MOs and upper, empty MOs, while towards d ~ 2 



practically all MOs are democratically occupied with half of the number of electrons they 
can accommodate. In the former case, electrons are completely delocalized over the whole 
nanoring, whereas in the latter, they are localized on QDs. Rephrasing in a more actual 
language, the quantum entanglement of electron states is nearly perfect in the limit of small 
interdot spacing (d > 1), while electrons become practically disentangled for large interdot 
spacing (d < 2). The curves for the von Neumann entropy S = — Trplog 2 p (p being the 
reduced density matrix) [47|, |48( , the quantity used to quantify quantum entanglement, pre- 
sented in Fig. [H nicely visualizes how the system smoothly evolves from the limit of perfect 
entanglement (S = 2) to that of perfect disentanglement (S = 1) by increasing d. For later 
purposes, it is also useful to monitor the weights p , pi, p 2 , ■ ■ ■ of the multielectronic config- 
urations in the wave function where none or several QDs are doubly occupied (X^P« = !)• 
The weights of these multielectronic configurations in the nonionized ground state $ are also 
shown in Fig. HJ Because electron hopping dominates at small d, one or even several QDs 
can be doubly occupied. On the contrary, on-site Coulomb repulsion U precludes double 
occupancy at larger d, where it is the dominant energy term, and all QDs tend to become 
singly occupied at half filling. 

Within the MO-approximation, Eq. ([2]) immediately yields wmo — 1 f° r au occupied 
electron Bloch states. Because Bloch states with +k and —k are degenerate for k ^ 0, N/2, 
out of the N/2 occupied Bloch states (for even N) there should be one contribution from the 
state with k = and (N — 2)/4 contributions from the states with \k\ ^ 0, i.e., (N + 2)/4 
distinct ionization signals within the MO. Therefore, to assess the validity of the single- 
particle picture, we shall (i) compare the number of spectral lines in the ionization spectrum 
with [N + 2)/4 and (ii) monitor the deviations from unity of the quantities w. 

IV. RESULTS 

In this section, we shall present in detail numerical results for ionization in half-filled 
nanorings consisting of six and ten QDs. Excepting for Fig. @fc, in all other cases of Sees. 
IIVAI and IIVBI we only depict the ionization signals with significant spectroscopic factors. 
In addition to these, there are numerous spectral lines of very small but still non-vanishing 
intensities, a fact to be discussed in Sec. fVTl 

Because of spatial symmetry, free (Bloch) electron states satisfy exactly the Hartree-Fock 
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FIG. 1: (Color online) Dependence on the interdot spacing d of the entropy (5) and the total 
weight po , pi, P2 etc. of the multielectronic configurations in the neutral ground state <3? of half- 
filled nanorings with no, one, two, etc. doubly occupied dots, respectively. For the nanoring with 
ten dots and electrons, the practically vanishing weight of the multielectronic configurations with 
five doubly occupied dots is not shown (p$ ~ 0). In both cases, for sufficiently distant QDs, the 
expansion of the neutral ground state is practically exhausted by multielectronic configurations for 
which none QD is doubly occupied (po ~ !)• The numbers of QDs in the nanoring are specified in 
the legend. 



(SCF) equations (see e.g., Ref. |49|). As a consequence, all ionization energies scale as the 
hopping integral to within the MO-approximation. Therefore, the curves for MO-ionization 
energies fall off exponentially with d. In the figures of this section, we shall not draw these 
rather trivial lines of the MO-approximation. We shall present instead the less trivial curves 
(depicted by black dashed lines in the panels showing the ionization energies) for the lowest 
ionization energy E\ oc in the large d limit, where electron hopping is negligible (to - ► 0) 
and none QD is doubly occupied. In this limit, the energy of the neutral ground state is 
Neh + NV . Removal of an electron from this state yields an ionized state with the energy 
(N — 1)eh + (N — 2)V. The difference of these energies represents the lowest ionization 
energy E\ oc = —Eh — 2V in the limit of perfect localization. 

Before going to discuss correlation effects on ionization spectra in detail, to give a flavor 
on the substantial difference between the MO and exact ionization spectra, we present an 
example in Fig. [2j Full spectra for all symmetries are shown in Fig. [2] for both six- and 
ten-dot nanorings, the two cases we shall later examine at length. There and in all of 
the other cases presented below, we have set Eh = — 4.504 eV, a choice that only fixes the 



8 



origin for energy [see Eq. (pQ)]. For a first glance assessment of the important role played by 
correlations for the spectra of Fig. [2] it suffices to say that, if the MO picture were valid, the 
spectra of Fig. [2h and [2b would only consist of two and three lines, respectively. 
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FIG. 2: (Color online) Exact ionization spectra for six- and ten-QD nanorings with interdot 
spacing d = 1.61. Notice that, in the MO spectra there exist only two and three lines (all with 
w = 1), respectively. The MO-lines are located at 4.556 eV and 4.609 eV in the former case, and 
at 4.536 eV, 4.589 eV, and 4.556 eV in the latter case. 



A. Six-dot nanorings 



Before discussing the numerical results for six-dot nanorings (point group D^), it is useful 
to note that, within the single-particle picture, there are six MOs. Ordered by increasing 
energy these are: a nondegenerate a\ g MO (k = 0), two degenerate e± u (k = ±1) and tig 
(k = ±2), and a nondegenerate b\ u MO (k — 3). At half filling the Oi g - and ei u -MOs are 
occupied by two and four electrons, respectively, whereas 6 lu - and e 29 -MOs are empty. 

From the multiplication table of the group Dg^ and because the neutral ground state 
possesses A± g symmetry, the ionization from an MO of a certain symmetry will bring the 
nanoring into a state of the same symmetry, e. g., removing an electron from an ei u -MO 
yields an ionized i?i u -state. Below, we shall discuss separately the ionization processes, 
depending on whether the MO wherefrom an electron is ejected would be occupied or not 
within the MO picture. 

Within the MO picture only two ionization processes are possible, because only the a l9 - 
and ei u -MOs are occupied. Accordingly, only two ionization processes (denoted by 1 in 



Figs. [3b, and [3fc, respectively) are significant at small d. The latter, which corresponds to 
the lowest ionization energy (HOMO-ionization), brings the nanoring into the ionized ground 



state. It is denoted by 1 in Fig. [3b and was already considered in Ref. [12|. With increasing 
d, the spectroscopic factors of these processes decrease from the MO-value w — 1. This 
represents one manifestation of electron correlations. Another manifestation is the fact that 
they become even significantly smaller than the corresponding MO occupancies: Wk,i < n^ 
(k = 0,1). (Henceforth we omitt the spin label whenever unncessary and write, e. g., rik 
instead of nk^-) Compare the curves denoted by 1 with the MO populations n = n^ l9 and 
rii = tie 1u in Figs. [3b and [3b, respectively. According to the sum rule (J3J), other eigenstates 
should aquire significant spectroscopic factors, and this is indeed what one observes in Figs. 
[3b and[3b. The differences 1— n lj0 can be considered a measure of the extent to which electron 
correlations affect the e\ u - and ai 9 -MOs of the neutral ground state. In comparison with the 
A l9 -ionization spectrum, where two ionized /l l9 -eigenstates practically exhaust all spectral 
weight in the whole <i-range of interest, that for Ei u -symmetry is somewhat richer; compare 
Fig. Eh with Fig. [3b- In the latter case we encounter an avoided crossing at d — 1.633, 
where the energy difference between the curves denoted by 4 and 3 attains its minimum 
values of 28.7 meV. Later on in this paper we shall return to a more detailed analysis of 
the phenomenon of avoided crossing. Here we only note that, by considering the diabatic 
approximation for the states participating to this avoided crossing, basically only three 
(diabatic) states give nonvanishing contributions to the E^-spectrum. 

Let us now switch to ionization from MOs that would be empty (in our case, e<ig and b\ u ) 
if the MO picture were valid. Such ionization processes are due to ground state correlations 



24j ; in the absence of the latter they do not obviously exist. Our results for these ionization 
processes are collected in Figs. H] and [5j In accord with the fact that electron correlations in 
the neutral ground state increase with d, a general trend of increasing spectroscopic factors 
is visible in these figures. Particularly interesting is the .E^-ionization spectrum, where a 
series of avoided crossings can be seen. As visible in Figs. 0b, one avoided crossing for each 
of the pairs of states (1, 2), (3, 4), and (5, 6) occurs at d ~ 1.765, d ~ 1.685, and d ~ 1.555, 
respectively. The corresponding minimum values of the energy differences amount 0.16 meV, 
1.0 meV, and 3.1 meV, respectively. In fact, in the avoided crossing at d ~ 1.555, not only 
two, but rather three (namely, 4, 5, and 6) states are involved. This fact is invisible in Fig.Hb, 
because the spectroscopic factor of the state 4 in this region is insignificant. In comparison 
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FIG. 3: (Color online) Ai g - and £a u -spectral factors w and ionization energies e versus interdot 
spacing d in six-QD nanorings. The numbers i = 1, 2, ... in the legend label the ionized eigenstates 
^k,i [cf. Eq. fl5J]. In panels (b) and (d), the black dashed line corresponds to the lowest ionization 
energy —en — 2V in the limit of perfect localization. 



with the i?2g~spectrum, the spectroscopic factors of the I?i u -spectrum are smaller: compare 
Fig. [5^ with UK This agrees with the fact that the lower the unoccupied MO, the more it 
becomes populated due to the ground state correlations. 

Actually, the states targeted by ionization in the case considered in this section represent 
eigenstates of a six-QD nanoring with five electrons. In this nanosystem, as discussed 
recently 14j], optical absorption allows to target eigenstates of A\ g , A2 g and E2 9 symmetries. 
So, ionization and optical absorption provide complementary information. Some states can 
be studied only by one method {E\ u and B\ u states only by ionization, A2 g states only by 
optical absorption), while other states (A\ g and E2 g ) by both methods. For the latter, the 
information extracted from optical absorption can be compared with that from ionization. 
By comparing the ionization and optical absorption spectra for A lff -symmetry (the present 



11 



Fig. [3^ and Fig. 6f of Ref. [l4|, noting the logarithmic scale in the latter), one can conclude 
that, in both cases, only the first two Ai 5 -states give an important contribution. Likewise, 
the same first six E^-states are the only ones that are important both for ionization and for 
optical absorption: compare the present Fig. [4^ with Fig. 6h of Ref. [l_4J. Still, all the three 
avoiding crossings at d ~ 1.765, d ~ 1.685, and d ~ 1.555 are visible in both spectra. The 
only difference is the participation in the optical absorption of three states at the avoided 
crossings at d ~ 1.685 and d ~ 1.555 (3, 4 and 5 in the former case, and 4, 5, and 6 in the 
latter), which is invisible in the E^-ionization spectrum. This is due to the insignificant 
spectroscopic factors of the states 3 and 4 in the corresponding regions (see Fig. Hh)- 

B. Ten-dot nanorings 

For rings with ten QDs, within the single-particle picture there are ten MOs. Ordered by 
increasing energy, they are: a nondegenerate a\ g MO (k = 0), two degenerate t\ u [k = ±1), 
C2 9 {k = ±2), e^u {k = ±3), e^ g {k = ±4) MOs, and a nondegenerate &i n MO {k = 5). 
Similar to the case analyzed in Sec. IIVA[ the ejection of an electron from an MO of a given 
symmetry leads to a final ionized state of the same symmetry, e. g., removing an electron 
from an e4 g -MO yields an ionized E 4s -state. This follows from the properties of the D 10i /j 
point group and the fact that the neutral ground state possesses A\ g symmetry. 

The numerical results for ionization in ten-dot nanorings are presented in Figs. [6] and H, 
separately for each symmetry. They confirm the general features already observed in six-QD 
nanorings. For the ionization processes having the symmetries (Ai g , E lu , and E 2g ) of the 
MOs completely occupied within the MO-picture (let us call them "occupied" MOs), there 
exist one signal (let us call it the "main" signal) per each symmetry whose spectroscopic 
factor is close to the ideal value wmo = 1 at d > 1. With increasing d, it decreases from 
wmo — 1 and becomes more and more smaller than the corresponding MO-population 
(n = riA lg , rii = n E lu , n 2 = n E 2g )- Concomitantly, a few more ("satellite") signals aquire 
more and more significant spectral weight, such that at larger d they can even dominate 
over the "main" signal. See Fig. [61 Interestingly, the deeper the "occupied" MO, the more 
is it affected by correlations at larger d; the relative weight of the "satellite" signals to 
the "main" signal for A^-symmetry (deepest "occupied" MO) is larger than that for E lM - 
symmetry, which is in turn larger than for E2 5 -symmetry (highest "occupied" molecular 
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FIG. 4: (Color online) -E^-spectral factors w (a) and ionization energies e (b) versus interdot 
spacing d in six-QD nanorings. The numbers i = 1, 2, ... in the legend label the ionized eigenstates 
^A:^ [c/. Eq. ([2"D]. In panel (b), the black dashed line corresponds to the lowest ionization energy 
—£h — 2V in the limit of perfect localization. Notice the numerous spectral lines with very small 
intensities in panel (c), which represents panel (a) redrawn using the logarithmic scale on the 
ordinate, demonstrating that the extended Hubbard model is characterized by a hidden quasi- 
symmetry (see Sec. I VII ) 



orbital, HOMO). 

At very small interdot distances {d > 1), the ionization from the "unoccupied" MOs (i. 
e., those unoccupied if the MO-picture were valid) is altogether ineffective. However, as 
a result of the strong electron correlations, the spectroscopic factors of these E 3u -, E± g -, 
and i?i u -ionization processes gradually increase with increasing d. Their spectral weight 
also becomes distributed over several lines. Similar to that from the "occupied" MOs, these 
spectroscopic factors are comparable among themselves, and also with those of the relevant 
ionization signals of the "occupied" MOs. See Figs. [7] and compare them with Figs. El 
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FIG. 5: (Color online) .Bi^-spectral factors w (a) and ionization energies e (b) versus interdot 
spacing d in six-QD nanorings. The numbers i = 1 and 2 in the legend label the ionized eigenstates 
^k,i [cf. Eq. (JS}]. In panel (b), the black dashed line corresponds to the lowest ionization energy 
—£h — 2V in the limit of perfect localization. 



To conclude, similar to the case of six-dot nanorings of Sec. IIVA| the MO-picture rapidly 
deteriorates with increasing d. 

From the inspection of the ionization factors of six-dot nanorings (Figs. [3], HI and EJ), one 
may be tempted to think that avoided crossings only occur for two-dimensional irreducible 
representations (Ei u and E 2g in that case) and not for the one-dimensional irreducible rep- 
resentations (Ai g and Bi u ). However, this is not true, as illustrated by the results for the 
one-dimensional irreducible representations of ten-dot nanorings presented in Figs. [6^ and 
[7^. For v4 ls -symmetry, two avoided crossings are visible in Figs. [6^ and b: one involving the 
states denoted by 3 and 4 at d ~ 1.561, another with the participation of the states 4 and 
5 at d ~ 1.670. The corresponding minimum energy differences are 5.33 meV and 0.26 meV, 
respectively. As seen in Figs. O and f, the avoided crossings for 5 lu -symmetry are more 
numerous. One of them, that at d ~ 1.643, even involves three states. 

Similar to the case of six-dot nanorings (Sec. II V Alb one can compare the information 
deduced by ionization for ionized nanorings with that obtained by optical absorption for 
nanorings with nine electrons and ten QDs. As discussed in Ref. [l4j, in the latter case one 
can target eigenstates of E\ u - and .E 3u -symmetry. 

Certain avoided crossings are visible both in ionization and in optical absorption spectra. 
This is the case for the avoided crossing at d ~ 1.755 involving the states 2 and 3, as well 
as the pretty broad one around d ~ 1.87. On the contrary, the dense avoided crossings 
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FIG. 6: (Color online) A\ g - } E\ u -, and E^^spectral factors w and ionization energies e versus 
interdot spacing d in ten-QD nanorings. The numbers i = 1, 2, ... in the legend label the ionized 
eigenstates ^k,i [cf- Eq. ([2"])]. In panel (d), the states 10 and 11 cannot be distinguished from the 
curves for 9 and 12 within the drawing accuracy. In panels (b), (d), and (f), the black dashed line 
corresponds to the lowest ionization energy —£h — 2V in the limit of perfect localization. 

occuring in the narrow range 1.723 < d < 1.737, which involve the states 9, 10, 11, and 12 
can be seen only in ionization (Fig. Eb and d), but not in optical absorption, because of the 
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FIG. 7: (Color online) E% u -, E^g-, and I?i u -spectral factors w and ionization energies e versus 
interdot spacing d in ten-QD nanorings. The numbers i = 1, 2, ... in the legend label the ionized 
eigenstates ^>k,i- The black dashed line corresponds to the lowest ionization energy — eh — 2V in 
the limit of perfect localization. 
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reduced intensities in the latter case; see Fig. 9d of Ref. 

The fact that more than two states can participate to an avoided crossing was already 
noted in the study on optical absorption [14]. An interesting situation can be seen in the 



1(3 



E^-spectrum (Fig. [TJi and b), where two avoided crossings are visible around d ~ 1.885: 
one at lower energies involves two states, another at higher energies with the participation 
of three states. These avoided crossings are i nsig nificant for optical absorption, because of 



the reduced spectral intensity of these states 14] 



A somewhat reversed situation occurs in the .E^-optical spectrum, where t 



re participa- 



tion of three states to the avoided crossing at d ~ 1.525 (see Fig. 9f of Ref. 1 141 ), while only 
two are significant (visible) in the E^-ionization spectrum of Fig. [7^. 

V. AVOIDED CROSSINGS 

The fact that avoided crossings involving energy curves of eigenstates of identical sym- 
metries represent a frequent phenomenon in metallic QD-nanorings has already noted in 
the study of optical absorption [14]. As a general characterization, as visible in Figs. Hk, Eh, 
Eb, Eb, [7h, [7b, and [7b they only occur at larger d, where correlations are important. Two (or 
more) states of identical symmetry described by a single Slater determinant, where the MOs 
are either occupied or empty, cannot come too close in energy; basically, their energy differ- 
ence is determined by the energy differences of their occupied MOs. Things change at larger 
d, where to ceases to be the dominant energy scale and frustration due to the U- and V-terms 
becomes important. For illustration, we present in Fig. [S] results for the pair of E 2g -states 
(E| and E\ ) of the six-QD nanoring (cf. Fig. HJ). As seen Fig. [8b, at the point d ~ 1.685 
the weights of the multielectronic configurations where none, one, or two QDs are doubly 
occupied rapidly interchange between the mates involved in the avoided crossing. The con- 
tributions to energy of the to-, U-, and V^-terms, i. e., — to ^2i a (^k,i\(a>icr a i+i o- + h- c -)l^fc,«)> 
U ^2i{$! k,i\ni tW, jj^j), and VJ2i('^k,i\ n i'ni+i\^k,i), respectively are sensitive to dot occupa- 
tions. Indeed, their curves at the avoided crossing, depicted in Fig. [Sb, behave accordingly. 

One may ask at this point whether avoided crossings only occur in the extended Hubbard 
ring, or also in the restricted (i. e., V = 0) Hubbard ring. The results presented in Fig. 
[91 which represents the counterpart of Fig. [8] for V = 0, reveal that they are also present 
in the latter case. Again, the physics behind avoided crossings is the interplay between the 
competing terms (in this case to- and [/-terms). For each eigenstate involved in an avoiding 
crossing, the individual (to- and £/"-) terms exhibit jumps at the avoided crossing point that 
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FIG. 8: (Color online) d-dependence of the weights Po,pi,P2 of the various multielectronic config- 
urations (a) and the separate contributions of the to~, U-, and V- terms to the total energy (b) 
of the ionized eigenstates E\ (thick lines) vs. E^ (thin lines). See the main text and the caption 
of Fig. [Q 

compensate each other in their sum (the total energy), which is represented by a smooth 
curve around this point. To avoid confusions, we note that the jumps at d ~ 1.495 visible 
in Fig. [9] for the curves of the state E\ are related to another avoided crossing between the 
states E\ and E\ . The latter, which is not shown in the figure, is characterized by jumps 
opposite to those for the i^f -curve. 




1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 d 




FIG. 9: (Color online) d-dependence of the weights Po,pi,P2 of the various multielectronic config- 
urations (a) and separate contributions of to - an d V- terms to the total energy (b) of the ionized 
eigenstates £f 9 (thick lines) vs. E^ (thin lines) in the case V = 0. See the main text and the 
caption of Fig. [Q 

From a pragmatic standpoint, one may wonder whether avoided crossings are of impor- 
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tance at all; since what one can measure there is merely a smooth bright diabatic state, and 
not two (or several) individual (adiabatic) states with identical symmetry exhibiting rapid 
variations. However, one should mention that this applies to the case of the ideal nanor- 
ing with equidistant QDs. In a distorted nanoring, there will be avoided crossings between 
the ionization energies of the different symmetries of the ideal nanoring. Besides, similar 
to ordinary molecules 3jJ, strong nonadiabatic effects can be expected in the presence of 
phonons. 



VI. HIDDEN QUASI SYMMETRY 



Most importantly, this investigation gives further support to the most intriguing aspect 
of electron correlations previously encountered in optical absorption 141 ]: in spite of the 
fact that correlations are strong, out of very numerous ionized eigenstates allowed by spatial 
symmetry to contribute to the ionized spectrum, only very few states possess a significant 
spectroscopic factor. 

The ionization spectra depicted in Figs. [2^ and b are indeed richer than those, which 
comprise two and three lines, respectively, obtained within the MO approach. However, 
the number of the ionization signals is much smaller than the number of states, which are 
expected to contribute in a strongly correlated system. To properly assess the intriguing 
aspect of correlations, one should compare the number of significant ionization signals with 
the substantially larger number of states allowed to contribute by spatial symmetry and spin 
conservation. For the six-dot nanoring of Figs. [2^, there are ten significant ionization signals 
but 178 allowed transitions. For the ten-dot nanorings, the number of 15 lines visible in 
Figs. [2b should be compared with the number ~ 10 4 of states allowed by spatial symmetry 
and spin conservation. 



In addition to well established symmetries 38 



39 



40 



41. 


42. 


43 



44 



45], the restricted 



(i. e., V = 0) Hubbard model for rings is known to possess hidden symmetries 461 ] . Similar 
to the case of optical absorption, the scarcity of the ionized spectra of nanorings described 
by the extended Hubbard model discussed here points towards a hidden guasi-symmetry: 
besides a few signals with significant spectroscopic factors, there are numerous lines of very 
small but definitely non-vanishing intensities. These almost vanishing signals, which exhibit 
a regular rf-dependence, can only be seen in Fig. lib, due to the logarithmic scale on the 
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ordinate. In all of the other figures presented so far, they are invisible within the drawing 
accuracy. Therefore, one may be tempted to think that this hidden quasi-symmetry will 
evolve into a true hidden symmetry as V — > 0. However, this is not the case, and for 
illustration we present in Fig. [TU] results for i? 2s -ionization in six-QD nanorings for V — 0, 
which represents the counterpart of Fig. H] at V 7^ 0. The presence in Fig. [TUb of numerous 
weak but non-vanishing spectral factors, also displaying a regular <i-dependence and being 
many magnitude orders larger than inherent numerical inaccuracies, demonstrates that the 
restricted Hubbard model is also characterized by a hidden g<iasz-symmetry, similar to the 
extended Hubbard model. 
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FIG. 10: (Color online) -E^g^spectral factors w (a) and ionization energies e (b) versus interdot 
spacing d in six-QD nanorings for V = (counterpart of Fig. H|). The numbers i = 1, 2, . . . in the 
legend label the ionized eigenstates *$>k,i [of- Eq. (|2j)]. Notice the numerous spectral lines with very 
small intensities in panel (c), which represents panel (a) redrawn using the logarithmic scale on 
the ordinate, demonstrating that the restricted Hubbard model is also characterized by a hidden 
guasi-symmetry. 
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Of course, the actual task of theory is to specify the hidden quasi-symmetry more pre- 
cisely. An interesting observation concerns the number of significant spectral lines per each 
symmetry. As already noted, the number iVj of significant ionization signals is larger than 
expected within the MO picture, or, rephrasing, larger than the number Nih of one-hole 
processes possible in the neutral self-consistent-field (SCF)-ground state. However, this 
number can be comprehended if one considers in addition to N\h the number N 2 h-ip of 
two-hole-one-particle processes possible in the neutral SCF-ground state, wherein the ex- 
cited electron occupies the MO just above the Fermi level. That is, the observation is that 
N, = N lh + N 2h _ lp . 

For instance, for a six-QD nanoring, there are two 1/i-processes, one with A\ g — and 
another with .E^-symmetry. The aforementioned 2/i-lp-processes are represented by the 
products E^ x £fj x E 2g , E^ x A{] x E 2g , A{] x E^ x E 2g , and A^ x A^ x E 2g . 
By using the multiplication rules of D 6/l , they lead to one A lg -, two E\ u -, four E 2g -, and 
two -Bi n -processes. In this way, we arrive at two A\ g -, three E\ u -, four E 2g -, and two B\ u - 
ionization processes. This is just what one observes in Figs. 3k,Eb,Hl and El respectively, if 
one considers the diabatic bright state at avoided crossings 501 ]. 

For ten-QD nanorings, there are three 1/i-processes, of Ai g -, Ei u -, and i^-symmetry. 
The analysis is similar but more tedious for the relevant 2/i-lp-processes. They are repre- 
sented by the products E^ 1 x E^ 1 x E 3u , E^ 1 x E^ 1 x E 3u , E 2g l x A^ x E 3u , E^ x E^ 1 x E 3u , 
E Tu x E \u x E 3u, E \u x A Tg x E 3u, A^ x E^ 1 x E 3u , A^ x E^ x E 3u , and A^ x A^ x E 3u . 
In view of the multiplication table of the point group Di ,h, one gets three Ai g -, five E± u -, 
five E 2g -, six E 3u -, six E 4g -, and three i?i„-processes altogether. Again, these numbers 
agree with the numbers of significant ionization signals visible in Figs. [6] and respectively. 

The fact that for six-QD nanorings we can compute exactly all eigenstates enables us 
to continuously vary the model parameters and find out the counterparts of the states with 
significant spectroscopic factors in the limit d > 1, where U, V -C 4£ and the MO-picture is 
reliable. The energies of the aforementioned lh- and 2 h-lp-pro cesses are given in Table Ufor 
the limit U, V -*• 0. For d -»• 1, one finds t Q ~ 1.49 eV , U/4t ~ 0.049, and V/4t ~ 0.036. 
By inspecting the curves for the ionization energies of Figs. [3b, [3H, HI and [5] one can see 
that they compare favorably to those of Tabled In this limit, the main effect of interaction 
is to split the degenerate energies by amounts of the order of U and V. This demonstrates 
that, indeed, the ionization processes with significant spectroscopic factors evolve into the 
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TABLE I: The energies of the lh— and 2/i-lp-processes described in the main text for six-QD 
nanorings in the limit U, V — ► 0. 



Symmetry 


Process 


Ionization Energy 


Degeneracy 


E\ u 


1/; 
2h-lp 


— £h + *o 
-eh + 4*o 


1 

2 


Alg 


lh 
2h-lp 


-E H + 2*0 
—£H + 3*0 


1 
1 


E 2 g 


2h-lp 
2h-lp 


— £ff + 3*o 

-e H + 5*o 


3 

1 


B\u 


2h-lp 


-en + 4*o 


2 



lh- and 2/i-lp-processes occurring in the uncorrelated case. 

For a quantitative description suggested by the above considerations, it appears most 
straightforwardly to employ Cfc i(J | $) andcfc )(T cj[ Cfc 2j(T3 |$) with \k\, \k 2 \ < kp and \ki\ = kp+1 
for constructing linear independent vectors with appropriate spin and spatial symmetries. 
Here, kp denotes the Fermi wave vector {kp = n for closed shell rings with N = 4n + 2, 



3- 



We have used them as 



n being an integer) and |$) the exact neutral ground state 
basis (sub) set for diagonalization, and found approximate eigenvectors ^k,i of the ionized 
nanoring. They have been utilized to compute an approximate ionization spectrum via Eq. 
(J2j). Apart from the fact that the number of lines is correctly obtained, this straightforward 
approach fails, however, to quantitatively reproduce the ionization spectra. 

In the subsequent attempt to reproduce the ionization spectra, we have generalized this 
approach, by lifting the above constraints imposed on k, fci, and k%. In a further effort, we 
have alternatively employed the exact ionized ground state \l/ and the dressed particle-hole 
excitations c k g Ck 2ya2 \^) to construct a basis subset for diagonalization. Both aforemen- 
tioned approximations are in the spirit of the treatment based on dressed particle-hole exci- 
tations developed in Ref. [14|, where it turned out to be an insightful approximate method to 
study optical absorption. Unfortunately, none of these two methods is able to quantitatively 
reproduce the exact ionization spectra of stro ngly correlated QD-nanorings. 

Similar to the case of optical absorption [Lj, ll4J, in spite of a variety of attempts of 
analyzing the numerical results, unfortunately, we cannot offer a quantitative explanation 
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of the hidden quasi-symmetry behind the scarcity of the ionization spectrum of strongly 
correlated nanorings. What we can do is only to tentatively speculate on the nature of this 
symmetry. Namely, based on the above considerations, we claim the existence of a one-to-one 
correspondence between the lh- and 2 h-lp-pio cesses possible in the SCF-neutral ground 
state and the number of significant lines in the ionization spectrum of strongly correlated 
half-filled QD-nano ring s described by the extended Hubbard model. The basic postulate of 



the Landau theory 52|, |53|, |54j is the one-to-one map between the low-energy excitations 
of noninteracting and interacting electron systems. Applied to ionization, Landau theory 
predicts a one-to-one correspondence between the ionization signals and the one-hole (lh) 
ionized states, a fact contradicted by the exact ionization spectra. Therefore, our hypothesis 
on the number of relevant lines represents an extension of Landau's basic idea. 

VII. CONCLUSION 

Because of their widely tunable properties, metallic QDs assembled in regular nanoarrays 
represent ideal controllable systems for bridging the regimes of weak and strong correlations. 
The present study of ionization in metallic QD nanorings confirms the important impact of 

il, UM, UJ$ • ft demonstrates 
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electron correlations found in a series of previous studies 12, 
that the MO-picture of ionization in tunable metallic QD-nanorings completely breaks down. 
This breakdown affects all MOs. In ordinary small molecules, it is possible to separate 
the ionization spectra of valence electrons in two distinct regions, related to outer- and 



inner-valence electrons 
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241 ] . In general, the MO-picture holds in the former but 
breaks down in the latter. Such a separation cannot be made in the metallic QD-nanorings 
investigated here, where strong correlations have impact on all of the MOs. As a limiting 
case thereof, even the HOMO-ionization is drastically affected by correlations. This fact, 
already pointed out in Ref. Il2j, contrasts to the case of ordinary molecules. Nevertheless, 
there exists a similarity between QD-nanorings and ordinary molecules: as discussed in Sec. 
IIVB| the higher "occupied" MOs are less affected by electron correlations than the lower 
ones. 

In ordinary molecules, weaker electron correlations manifest themselves in ionization 



spectra as satellite lines of small intensities accompanying the main ionization signals 
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20 



y 



13, 



2l| . The main, more intense lines are the result of one-hole (lh) processes related 
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to the ejection of an electron from an MO. The satellite, less intense lines are related to 



excitations accompanying the main ionization . 



, often consisting of two-hole-one-particle 



(2h-lp) processes. In the case of strong correlations in molecules the intensity is distributed 
over numerous lines with comparable inten sity . This effect was termed the breakdown of 



the molecular orbital picture of ionization 22|, |23j, |24| . Main lines cannot be identified in 



this case, because the properties of an ordinary molecule are not tunable, and tracing back 
to the uncorrelated limit is impossible. 

The distinctive feature found in metallic QD-nanorings studied here is that, as visible in 
Figs. [3] and EJ in the presence of strong correlations the signals originating from the main lines 
of the uncorrelated case progressively lose intensity and become dominated by those that 
trace back to 2h-lp processes in uncorrelated nanorings. As concerns the number of lines 
in the ionization spectra, it remains astonishin gly scarce. Along with the similar results 
of the recent studies on optical absorption [13, ll4(, this finding gives further support to 
the existence of a hidden quasi-symmetry in the metallic strongly correlated QD-nanorings 
described within the extended Hubbard model. As shown here, it remains a quasi-symmetry 
even in the case of the restricted Hubbard model (V = 0). 

With regards to the hidden quasi-symmetry, amply documented here and in other recent 



works 
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14j by detailed numerical results, but for which a physical explanation is not yet 
available, we make the following remark. Exact numerical results played an important role in 
many cases, even including the Hubbard model. The fact that the one-dimensional restricted 
Hubbard model possesses accidental degenerqacies was amply discussed in the literature in 
the past. What initiated this issue were the exact numerical results published long time ago 



for benzene-like rings [56(. Although a physical explanation could not be given that time, 
the finding on accidental degeneracies became part of wisdom for the Hubbard community. 
It was only much later that a demonstration that these "accidental" degeneracies follow 



from certain nontrivial conservation laws could be given [57j. Similarly, we hope that our 
numerical findings will stimulate further investigations to unravel the nature of the hidden 
quasi-symmetry put forward here. 

We think that, the results of the present theoretical study on ionization along with those 
on optical absorption [13, [l4| are of interest and hope that they will stimulate scientists to 
fabricate and investigate metallic QD-nanorings. 



24 



Acknowledgements 

The authors acknokledge with thanks the financial support provided by the Deutsche 
Forschungsgemeinschaft (DFG). 



[1] J. Heath, C. Knobler, and D. Leff, J. Phys. Chem. B 101, 189 (1997). 

[2] C. P. Collier, R. J. Saykally, J. J. Shiang, S. E. Henrichs, and J. R. Heath, Science 277, 1978 

(1997). 
[3] R. Ingram et al., J. Amer. Chem. Soc. 119, 9279 (1997). 

[4] G. Markovich, C. P. Collier, and J. R. Heath, Phys. Rev. Lett. 80, 3807 (1998). 
[5] J. Shiang, J. Heath, C. Collier, and R. Saykally, J. Phys. Chem. B 102, 3425 (1998). 
[6] S. Chen et al., Science 280, 2098 (1998). 
[7] G. Medeiros-Ribeiro, D. A. A. Ohlberg, R. S. Williams, and J. R. Heath, Phys. Rev. B 59, 

1633 (1999). 
[8] S. Henrichs, C. Collier, R. Saykally, Y. Shen, and J. Heath, J. Amer. Chem. Soc. 122, 4077 

(2000). 
[9] J. Sampaio, K. Beverly, and J. Heath, J. Phys. Chem. B 105, 8797 (2001). 
[10] K. Beverly, J. Sampaio, and J. Heath, J. Phys. Chem. B 106, 2131 (2002). 
[11] M. A. Kastner, Physics Today 46, 24 (1993). 

[12] I. Baldea and L. S. Cederbaum, Phys. Rev. Lett. 89, 133003 (2002). 
[13] I. Baldea and L. S. Cederbaum, Phys. Rev. B 75, 125323 (2007). 
[14] I. Baldea and L. S. Cederbaum, Phys. Rev. B 77, 165339 (2008). 

[15] I. Baldea, A. K. Gupta, L. S. Cederbaum, and N. Moiseyev, Phys. Rev. B 69, 245311 (2004). 
[16] K. Siegbahn et al., ESC A- Atomic Molecular and Solid State Structures Studies by Means of 

Electron Spectroscopy, North-Holland, Amsterdam, 1967. 
[17] K. Siegbahn et al., ESCA-Applied to Free Molecules, North-Holland, Amsterdam, 1969. 
[18] K. Siegbahn and L. Karlsson, in Handbuch der Physik, edited by S. Fliigge, page 215, Springer, 

Heidelberg, 1982. 
[19] D. W. Turner, C. Baker, A. D. Baker, and C. R. Brundle, Molecular Photoelectron Spec- 
troscopy, Wiley, New York, 1970. 



25 



[20] J. W. Rabalais, Principles of Ultraviolet Photoelectron Spectroscopy, Wiley, New York, 1970. 

[21] J. Berkowitz, Photoab sorption, Photoionization and Photoelectron Spectroscopy, Academic 
Press, New York, 1979. 

[22] L. S. Cederbaum, J. Schirmer, W. Domcke, and W. von Niessen, J. Phys. B 10, L549 (1977). 

[23] L. S. Cederbaum, W. Domcke, J. Schirmer, and W. von Niessen, Physica Scripta 21, 481 
(1980). 

[24] L. S. Cederbaum, W. Domcke, J. Schirmer, and W. von Niessen, Adv. Chem. Phys. 65, 115 
(1986). 

[25] F. Remade and R. D. Levine, J. Amer. Chem. Soc. 122, 4048 (2000). 

[26] F. Remade et al., J. Phys. Chem. B 102, 7727 (1998). 

[27] F. Remade, C. P. Collier, J. R. Heath, and R. D. Levine, Chem. Phys. Lett. 291, 453 (1998). 

[28] F. Remade and R. Levine, J. Phys. Chem. B 105, 2153 (2001). 

[29] F. Remade, K. Beverly, J. Heath, and R. Levine, J. Phys. Chem. B 106, 4116 (2002). 

[30] One should remark here that not all the symmetries of the group Djvh are spanned by the 
MO-symmetries. This is the case, for example, of the eigenstates with A<i g or P>2 U of a six-dot 
nanoring, which cannot therefore be studied by photoionization. 

[31] H. Koppel, W. Domcke, and L. S. Cederbaum, Adv. Chem. Phys. 57, 59 (1984). 

[32] I. Baldea, H. Koppel, and L. S. Cederbaum, Phys. Rev. B 55, 1481 (1997). 

[33] I. Baldea, H. Koppel, and L. S. Cederbaum, Phys. Rev. B 69, 075307 (2004). 

[34] D. Bullet, R. Haydock, V. Heine, and M. J. Kelly, in Solid State Physics, edited by H. Erhen- 
reich, F. Seitz, and D. Turnbull, Academic, New York, 1980. 

[35] P. Fulde, Electron correlations in molecules and solids, in Springer Series in Solid-State 
Sciences, volume 100, Springer- Verlag (Berlin, Heidelberg, New York), 1991. 

[36] E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

[37] The fact that \^k,a) = Ck,a\&) is a state of well defined S z {= ^1/2) and S(= 1/2) straight- 
forwardly follows from the commutation relations of Ck )0 - an d the spin operators, and the fact 
that S Z \Q) = 5+| ( ! > ) = (<£ is a singlet). The commutation relation [-S^c^jJ = 2 c fc,| yields 
Sz\^k,i) = ~^~h\^k,i)- Similarly, because c/%^ and the spin raising operator S+ commute, one 
immediately gets S' + |* fc ,|) = ° and thence 5 2 |* fe; |) = {S 2 Z + S-S+ + S z )\^ k ^) = ||^fc,|). 

[38] A. Natan, Integrable models in condensed matter physics, in Series on Modern Condensed 
Matter Physics, Lecture Notes of ICTP Summer Course, edited by S. Lundquist, G. Morandi, 

26 



and Y. Lu, volume 6, page 458, World Scientific Singapore, 1992. 



[39] A. Natan, Integrable models in condensed matter physics, 1994, lcond-mat/9408101 



[40] E. H. Lieb and F. Y. Wu, Physica A: Statistical Mechanics and its Applications 321, 1 (2003). 

[41] B. S. Shastry, Phys. Rev. Lett. 56, 1529 (1986). 

[42] B. Shastry, J. Stat. Phys. 50, 57 (1988). 

[43] H. Grosse, Lett. Math. Phys. 18, 151 (1989). 

[44] M. Wadati, E. Olmedilla, and Y. Akutsu, J. Phys. Soc. Jpn. 56, 1340 (1987). 

[45] E. Olmedilla and M. Wadati, Phys. Rev. Lett. 60, 1595 (1988). 

[46] P. B. Ramos and M. J. Martins, J. Phys. A: Mathematical and General 30, L195 (1997). 

[47] C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schumacher, Phys. Rev. A 53, 2046 (1996). 

[48] G. C. Ghirardi and L. Marinatto, Phys. Rev. A 70, 012109 (2004). 

[49] A. L. Fetter and J. D. Walecka, Quantum Theory of Many Particle Systems, McGraw Hill, 

New York, 1971. 
[50] The bright diabatic state approximates the mate with significant spectroscopic factor of a pair 

of states involved in an avoided crossing. For instance, in Fig. [H it describes the state 1 for 

d < 1.765 and the state 2 for d > 1.765. 
[51] Unless d is very close to unity, using the SCF-ground state \&scf) instead of |<3?) is a very 

poor approximation, a fact easily understandable in view of the strong correlations. 
[52] L. D. Landau, Sov. Phys. JETP 3, 920 (1956). 
[53] L. D. Landau, Sov. Phys. JETP 5, 101 (1957). 
[54] L. D. Landau, Sov. Phys. JETP 8, 70 (1959). 
[55] T. Aberg, Phys. Rev. A 2, 1726 (1970). 

[56] O. J. Heilmann and E. H. Lieb, Annals of the New York Academy of Sciences 172, 583 (1971). 
[57] E. A. Yuzbashyan, B. L. Altshuler, and B. S. Shastry, J. Phys. A: Mathematical and General 

35, 7525 (2002). 



27 



